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A quantum mechanical analysis of the bremsstrahlung in a decay of 210 Po is performed in close 
reference to a semiclassical theory. We clarify the contribution from the tunneling, mixed, outside 
barrier regions and from the wall of the inner potential well to the final spectral distribution, 
and discuss their interplay. We also comment on the validity of semiclassical calculations, and 
the possibility to eliminate the ambiguity in the nuclear potential between the alpha particle and 
daughter nucleus using the bremsstrahlung spectrum. 
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■ Bremsstrahlung during a decay and fission is one of the intriguing current topics of nuclear physics. Kasagi et 
ON ■ al. measured the bremsstrahlung in the a decay from Po, and showed that the emission of low energy photons is 
systematically hindered from what one expects from the pure Coulomb acceleration of the a particle in the outside 
Oh' region of the potential barrier. They speculated that this reduction is caused by a destructive interference between 
the radiations in the tunneling and classically allowed regions, and thus might offer a possibility to learn about the 
tunneling time. They also calculated the spectral distribution based on a semiclassical theory of Dyakonov and Gornyi 
Q . Their calculations reproduce the experimental cross section well for low energy photons. In addition, an interesting 
thing is that their calculations give a hump in agreement with the measurement, which however has large error bars 
in that region. On the other hand, the intensity of the radiation observed by D'Arrigo et al. || for the alpha decay 
from 226 Ra and 214 Po is larger than the prediction of the pure Coulomb acceleration model. Also, their data show no 
structure in the spectral distribution as expected from a classical formula. In nuclear fission, a measurement by Luke 
, et al. gave an upper bound to the bremsstrahlung rate for the spontaneous fission of 252 Cf. 
^\ ' Recently, Papenbrock and Bertsch performed a quantum mechanical calculation for the bremsstrahlung in a decay 
in perturbation theory Their spectral distribution is monotonic, though it is within the error bars of the data 
OO ■ in Ref. Q. The authors claim that the contribution from the tunneling wave function under the barrier is small. 
ON ■ However, the definition of either classical or tunneling is not unique. More studies from various points of view are 
i clearly needed in order to clarify the situation. The aim of this paper is to resolve these puzzling problems relating to 
the spectral distribution, especially on the role of the tunneling region. To this end, we perform a quantum mechanical 
analysis in a way which has a close relation to semiclassical theory. As we show later, quantum mechanical calculations 
are required to have quantitative accuracy, while a semiclassical theory provides a clear physical understanding of the 
phenomena by bridging the quantum mechanical and classical calculations of the bremsstrahlung. It also provides 
a more definitive understanding of the role of quantum tunneling by naturally dividing the whole region into the 
classical, mixed and tunneling regions through classical turning points. 

Fermi's golden rule gives the following expression for the photon emission per photon energy during a decay in the 
dipole approximation, 

where V(r) is the interaction between the a particle and the daughter nucleus as a function of their separation 
distance r, m is their reduced mass and E 1 is the photon energy. The effective charge Z e g is given by Z e g — 
(Ai)Z a —A a Zu)/Ap, where Ad, Zd, A a , Z a and Ap are the mass and atomic numbers of the daughter nucleus and of 
a particle, and the mass number of the parent nucleus, respectively. <5>i(r) and are the appropriately normalized 

radial wave functions of the initial and final states of the a particle, respectively H. 

As in Ref. 0, we take the following simple model for the interaction between the a particle and the daughter 
nucleus 

V(r) = {Z a Z D e 2 /r)e{r - r„) - V e(r - r). (2) 
The initial and final state wave functions are given by 
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where Ri(r) — Go(r)i,kir) + iFo(rji,kir) and Rf(r) — Fi(rjf, kfr) cosa + Gi(r/f, kfr) sina, a being the phase shift. 
The wave numbers in these equations are given by, 



K l ^fi- 1 ^2m{E l + V a )- h = h- 1 ^/2mE l (5) 



K f ^hr x ^2m{E s + V Q ); k f = hr 1 ^2^E~ s (6) 

where E{ = E a , E a being the energy of the a decay, Et = E a — E~ = E a — Tuo and rj the Sommcrfeld parameter. 
The transition matrix in Eq.(Q) now reads, 
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where A(r) stands for d r V = — Z a Z]je 2 r^ 2 , and the wall contribution lyj is given by 

I w = (Z a Z D e 2 /r Q + V )R f (r Q )R l (r ) (8) 

If we denote the external turning points in the initial and final states by r e i and r e / , respectively, the integration 
in Eq.(Q) can be naturally divided into three parts, i.e. the integrations between ro and r e i, between r e i and r e f 
and between r e f and oo. We call them the tunneling, mixed and outside regions, respectively. Straight-forward 
calculations of each integral using the exact Coulomb wave functions will clarify the role of each term, especially 
the role of the tunneling region. Before we present the results of such quantum mechanical calculations, we wish to 
present semiclassical formulae for radiation. 

They can be derived by replacing the quantum mechanical Coulomb wave functions by their semiclassical represen- 
tations, which read [|| 
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where r e is the external classical turning point, and ki and Kg are the wave numbers for the orbital angular momentum 
I in the classically allowed and forbidden regions, respectively. Notice that Gg and Eg increases and decreases 
exponentially, respectively, as r moves deeper inside the tunneling region. We remark that the phase shift a in Eq.(^) 
is proportional to the tunneling probability and is very small. This can be proved by matching the semiclassical wave 
function under the potential barrier to the wave function in the region of the potential well. We thus express tan a as 

tana = C a P t f, P t f = e Jr o , (11) 



where P t f is the tunneling probability in the final state in the WKB approximation. In Eq.Ql 1|) and in what follows, 
the meaning of the lower indices will be obvious, e.g. K\f is the wave number with angular momentum 1 in the final 
state. 

We now derive the semiclassical formulae for the contribution to the integral in Eq.(^) from the tunneling, mixed 
and outside regions separately. 

1. Tunneling region: The main contribution from the tunneling region comes from the F\ -Gq term, whose semiclassical 
representation is given by 
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In obtaining the last expression, we expanded the Ki/(r) — Koi( r ) U P to the leading order of ^7f7~7jp + ^ and 
introduced the time parameter along the (imaginary) time axis r% by 
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li = Ti(r e i) is the tunneling time in the initial channel. The lower index i of r, shows that the time is related to the 
distance r through the velocity in the initial state. Notice that the radiation amplitude in this region is given by the 
Laplace transform of the acceleration if one could ignore the centrifugal term. The G\ ■ Go term is the largest as long 
as the integral itself is concerned. However, the contribution of this term is reduced by the small value of tana in 
front of G\. We have a similar expression for the G\ ■ F term, which is also expected to be small. The Fx ■ F term 
is expected to be the smallest. 

2. Mixed region: The mixed region, which is classically allowed in the initial state but forbidden after radiation, is 
beyond the scope of classical calculations in Refs. jl| and p|. Although quantum mechanical calculations are required 
to quantitatively estimate the contribution from this region as we show later, it is still interesting to see its semiclassical 
representation. The contribution of the Fx term which dominates in this region reads, 

f' S A(r) Fl ( Vf ,k f r) Rl (r)dr ~ \ f" A(r) ( - ffi ) e " dr (l 5 ) 

We show later that this region strongly influences the behaviour of the spectral distribution at high energies. 

3. Outside the barrier: The main contribution in the classical region, i.e. outside the barrier, comes from the Fx term. 
Its semiclassical expression reads, 

i t A{r) Fl{ ^ k * r) Rdr)dr ~ ~^r e " J A ^ )] \k^m ) e f (16) 

This is nothing but the classical formula for the bremsstrahlung except for the existence of the centrifugal potential 
term. The lower index / of tf means that the time t is related to the distance r through the velocity in the final 
state. We have ignored the term where the action integrals in the initial and final states enter with the same sign, 
because the sign of the integrand in that case rapidly changes as a function of r and the value of integration becomes 
negligibly small. 

We next consider the contribution from the wall. This term has been overlooked in the calculations reported in Ref. 
. As we clearly demonstrate later, this is one of the key issues to resolve the discrepancy between the calculations 
in Refs. jlj] and 0. The contribution of the Fq term can be safely ignored unless the energy of the a particle is very 
close to the top of the potential barrier. Using Eq. (pr|) , the remaining terms can be expressed as, 

Rf{r )G Q ( m , h r«)~\{ ( \ h -, . )La(l + 2C B )e-C^V*Jr^^, (17) 
2\Kif(r )Koi(r )J 

We now analyse the bremsstrahlung in the a decay from 210 Po, where E a =5.3 MeV. We determine the potential 
parameters in the same way as in Ref. ||, i.e. by matching the quantum mechanical wave functions at ro, and by 
calculating the decay width by normalizing the current at r e to the probability inside. An alternative is to use the 
semiclassical quantization rule and the R-matrix theory for the decay rate (see Eqs.(3) and (4) in Ref. 0). We 
found, however, that the semi-classical quantization rule is not valid with significant influences on the bremsstrahlung 
spectrum (see later). 

Figure 1 shows the spectral distribution calculated quantum mechanically with one of the potential parameter sets, 
Vq = 21.37 MeV and ro = 8.055 fm, which gives 5 nodes to the wave function for the radial motion. The general 
behaviour of the theoretical results do not depend so much on the particular choice of the potential parameters (see 
below, however). This potential parameter set corresponds to that in Refs. (|J^]. The thick solid line represents the 
total photon emission probability. It is a monotonically decreasing function and is consistent with the data except 
for the last data point. The figure also shows that various components have nearly the same order of magnitude, and 
that the final spectrum is the result of a complicated interference among them. In order to see the role of tunneling 
and classical regions and the interplay of each contribution more clearly, we show in Fig. 2 specral distributions by 
dividing into the classical and tunneling contributions. We present two different combinations depending on whether 
we consider the mixed region classical or tunneling. We call the part containing the outside barrier component the 
classical. In either case, the soft photon emission is dominated by classical contributions. This follows from the 
strong cancellation between the wall contribution and the contribution from the integral in the tunneling region. 
If the potential has a finite slope at the inner turning point and if one could ignore the centrifugal term, the net 
contribution from the tunneling region should be zero in the soft photon limit. Our result matches with this general 
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aspect. For high energy photons, the classical and the tunneling contributions have nearly the same magnitude and 
interfere destructively leading to a much smaller photon emission than the the prediction of the classical theory. 
The dot-dashed line in Fig. 2 shows the photon emission probability obtained by summing the contributions from the 
outside barrier and tunneling regions. Interestingly, this has a hump at high photon energies, just like that discussed 
in Ref. [jD, which overlooked the wall contribution and ignored the contribution from the mixed region. However, this 
interesting interference pattern is hidden by these contributions. 

As well known, the resonance position and the decay width do not uniquely determine the potential. Each po- 
tential gives a different value for the coefficient C a in our formalism. An interesting question is then whether the 
bremsstrahlung can be used to give additional constraints to the potential. If one restricts to the simple square well 
potential in Eq.(||), various sets correspond to different number of nodes n in the radial wave function inside the 
potential well. We found that the spectral distribution is almost indistinguishable as long as n < 7. For larger values 
of n, one starts to notice a difference in that the soft photon emission probability is predicted to be recognizably 
larger than the prediction of the classical theory. Note that the soft photons mean here those whose energy is, say, 
smaller than 100 keV. The emission probability of zero energy photons is not influenced by the potential parameters 
as discussed in Ref. ||. As an example, we show in Fig. 3 the spectral distribution calculated by the potential set, 
Vq = 107.3 MeV and tq — 7.947 fm, which corresponds to n = 11 and matches with a simple estimate based on the 
Pauli principle. In this respect, it would be extremely interesting to extend the data towards lower photon energies 
to examine whether the shallow or the deep potential agrees better with the data of bremsstrahlung. C' a gradually 
decreases from 0.041 to —0.055, and from —0.25 to —0.28 for the potentials in Figs.l and 3, respectively, as E 1 changes 
from to 600 keV. 

We finally comment on the validity of semiclassical calculations. We already mentioned a problem with the semiclas- 
sical quantization rule in determining the potential parameters. Since C a is sensitive to the slope of the wave function 
at the matching point, it causes a serious error in properly describing the bremsstarhlung spectrum. We found that 
the wall contribution itself gets very small because of the cancellation between the F\ and G\ terms in Eq.([l7|) if we 
adopt any potential set determined by using the semiclassical quantization rule. Consequently, the contribution from 
the tunneling region dominates the spectrum at soft photons giving larger cross section than the classical prediction. 
Another problem is that two turning points r e i and r e / lie too close to each other to use naive semiclassical wave func- 
tions. The crosses in Fig.l represent the contribution from the mixed region calculated semiclassically. The deviation 
from the dotted line clearly shows the failure of the semiclassical calculations. Since the contribution from the mixed 
region strongly influences high energy spectrum, this is a serious problem. The semiclassical calculations reproduce 
the qualitative behaviour of the quantum mechanical results quite well for the contributions from the regions under 
and outside the barrier, but have some quantitative inaccuracy especially at high energies. 

In summary, our analysis shows that the final bremsstrahlung spectrum results from a subtle interferences of the 
contributions of the tunneling, mixed and classical regions as well as the wall of the the potential well, each of which 
has a comparable magnitude. Semiclassical as well as classical theories seem not to be reliable for describing these 
subtle interference effects, though they give some clear understanding of the phenomena. It will be very interesting if 
one could perform more exclusive experiments to pick up each contribution separately. The extention of data to more 
soft photon side is awaited to provide a stringent test of the potential between the alpha particle and the daughter 
nucleus. Our study is based on the assumption of the validity of the potential model. Though this assumption 
seems to work in our analysis, it is a very interesting question to examine more in detail whether the bremsstrahlung 
spectrum contains some information beyond the scope of the potential model such as the validity of the R-matrix 
theory and the preformation factor. 
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Figure Captions: 

Fig. l:Total photon emission probability and its decomposition. The potential parameters are Vq = 21.37 MeV and 
ro = 8.055 fm. The crosses are the mixed region contribution calculated semiclassically. All the other lines have been 
obtained quantum mechanically. The data are from Ref. 0. 

Fig. 2: Grouping of each component into classical and tunneling contributions. The potential parameters are the 
same as those for Fig. 1. 

Fig. 3:The same as Fig.2, but for V = 107.3 MeV and r = 7.947 fm. 
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